Pseudogap and spin fluctuations in the normal state of electron-doped cuprates 
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We present reliable many-body calculations for the t-t'-t"-U Hubbard model that explain in detail 
the results of recent angle-resolved photoemission experiments on electron-doped high-temperature 
superconductors. The origin of the pseudogap is traced to two-dimensional antiferromagnetic spin 
fluctuations whose calculated temperature dependent correlation length also agrees with recent 
neutron scattering measurements. We make specific predictions for photoemission, for neutron 
scattering and for the phase diagram. 
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High-temperature superconductors (HTSC) still 
present one of the main contemporary challenges to 
condensed matter physics. While single-particle exci- 
tations, in the standard theory, are described by the 
quasiparticle (Fermi liquid) concept, that approach fails 
in HTSC. Quasiparticles can disappear altogether on 
certain segments of the would-be Fermi surface p| . This 
in turn leads to anomalous properties in transport and 
in thermodynamic data that are collectively referred 
to as pseudogap phenomena. Pseudogap phenomena 
occur both in hole-doped (h-doped) cuprates and in the 
electron-doped (e-doped) ones. Although the majority 
of research has focused on h-doped materials, experi- 
mental activity on e-doped cuprates has been steadily 
increasing in the last few years, as the community 
realizes the importance of looking at the big picture 
for HTSC. Strong electron-electron interactions and 
low dimensionality are the main stumbling blocks for 
theories of HTSC. In this context, one can argue that 
it is the lack of controlled theoretical approximations 
and consequent lack of quantitative predictions and 
detailed agreement with experiment that are the main 
reasons why there is no consensus on the correct theory 
in this field. We will show in this paper that the 
somewhat weaker coupling in e-doped HTSC leads to an 
unprecedented opportunity to obtain detailed agreement 
between theory and experiment. 

We use the single-band Hubbard model on a square 
lattice, which contains a repulsive local interaction U 
and a kinetic contribution of bandwidth W that is fitted 
to band structure using nearest t, second-nearest t', and 
third-nearest t" neighbor hoppings. Reliable calculations 
for the Hubbard model can be done when the expansion 
parameter U/W becomes less than unity. As was pointed 
out by several groups 0[|[3 Q , this is precisely the situ- 
ation that occurs in e-doped HTSC as doping is increased 
towards optimal doping. Smaller values of U/W physi- 
cally come from better screening in e-doped systems [5|- 
Note that the pseudogap mechanism discussed here can 
be different from that occuring at strong coupling 

We calculate the momentum k and energy lu depen- 



dent single-particle spectral weight A(k,Lo) and find that 
it is in detailed agreement with Angle Resolved Photo- 
Emission Spectroscopy (ARPES) on Nd 2 _xCe 2; Cu04 Q, 
EJQ- In particular, we explain the pseudogap observed 
in e-doped cuprates at binding energies of the order of 
a few hundred meV. To show that we correctly handle 
the physics of antiferromagnetic (AFM) fluctuations that 
is behind that pseudogap, we successfully compare our 
temperature-dependent correlation length to that mea- 
sured by neutron scattering on Nd2- a; Ce a ;Cu04 H, 
None of the recent theoretical works on e-doped cuprates 
S [J HI EJ El have obtained momentum, frequency 
as well as temperature dependent properties that are in 
detailed agreement with both ARPES and neutron data. 
We also find good estimates for the Neel temperature @ 
and for the d-wave superconducting transition temper- 
ature 12j . A few of our predictions are that: (a) The 
ARPES pseudogap disappears with increasing tempera- 
ture when the AFM correlation length £ becomes smaller 
than the single-particle thermal de Broglie wavelength 
£th. The corresponding pseudogap temperature T* is 
close to that found in optical experiments [l^. (b) In 
the pseudogap regime, the characteristic spin fluctuation 
energy in neutron scattering experiments is smaller than 
the thermal energy (renormalized classical regime) and 
the spin fluctuations are overdamped near the pseudogap 
temperature, (c) There is no large energy (> 100 meV) 
pseudogap at dopings larger than that of the quantum 
critical point (QCP), whose location we also calculate 
and which coincides with that observed experimentally 
for the antiferromagnetic to paramagnetic QCP ja, ^| . 

Since we are in a non-perturbative regime, even with 
U/W < 1, calculations are still a challenge. We use the 
non-perturbative Two-Particle Self-Consistent approach 
(TPSC) valid in that regime. The reasons are as fol- 
lows. TPSC satisfies at the same time the Pauli princi- 
ple, the Mermin- Wagner theorem and has a self-energy 
that is consistent with the sum-rule relating potential 
energy to single-particle properties 0, 0, 0, 0] . By 
contrast with other approaches Q| , TPSC also automati- 
cally takes into account quantum fluctuations that renor- 
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FIG. 1: Energy distribution curves A < (k,ui) = A(k,u)f(u) 
along the Fermi surface shown in the insets for (a) n = 1.15, 
U = 5.75* (b) n = 1.10, U = 6.25*. Lines are shifted by 
a constant for clarity. Spectra along two other directions for 
n = 1.15 are also shown: (c) (0, 0) — (n, n) and (d) (0.657T, 0) — 
(0.65-7T, 7r). See also the corresponding experimental data in 
Fig. 2 of Ref. Q and Fig. 1 of Ref. |(|. Band parameters are 
t' = -0.175*, t" = 0.05*, * = 350 meV while T = t/20. 



malize U (Kanamori screening). Results for the charge 
and spin structure factors, spin susceptibility and dou- 
ble occupancy obtained with the TPSC scheme satisfy 
conservation laws. As a theory then, TPSC fulfills cru- 
cial physical constraints. In addition, the above quan- 
tities have been found to be in quantitative agreement 
with Quantum Monte Carlo (QMC) and other numeri- 
cal methods [TR 111 \rA 111 ITal^ for both the nearest- 
neighbor [l5tll6| and next-nearest-neighbor |0| Hubbard 
models in two dimensions. In particular, the agreement 
with QMC for single-particle properties has been veri- 
fied 0, 0, |2(j even as we enter the pseudogap regime 
at low enough temperature. An internal accuracy check 
(l(| allows one to decide on the validity of the results in 
cases where QMC or other exact results are not avail- 
able as references. We do not consider the case of very 
large correlation length, where the quantitative details 
of TPSC become less accurate. Therefore, we are always 
within the domain where the approach is reliable. 

a. Energy dispersion curves One of the most strik- 
ing features of ARPES results in optimally e-doped 
cuprates is the disappearance of u> = excitations near 
the intersection of the non-interacting Fermi surface with 
the AFM zone boundary. These pseudogaped points in 
the Brillouin zone are called hot spots. Band struc- 
ture parameters that are input in the calculation lead 
to intersection of the Fermi surface with the AFM zone 
boundary. Throughout the paper, we present results 
for *' = —0.175* and t" = 0.05*, but very similar re- 
sults are obtained for another choice of parameters 3 
t' = -0.275* and t" = 0.0 |2l| . There are no phenomeno- 



logical parameters, such as mode coupling constants or 
additional renormalizations of U, as can occur in other 
theories. For actual comparisons with experiment, we 
take t — 350 meV since this allows us to find agree- 
ment with the experimentally observed overall ARPES 
dispersion. 

We present in Fig. ^ the single-particle spectral weight 
for two fillings, n = 1.15 and 1.10. Here, A(k,oj) 
is multiplied by the Fermi-Dirac distribution function 
/(lu) to allow direct comparison with experiment on 
Nd2_ a: Ce :c Cu04 (n= l + i) 1| One should also 
recall that experimental data contain a background at 
large binding energy whose origin is most probably ex- 
trinsic 23]. The results that we present in the next two 
figures, Figs. n an d 121 a- b)> are at temperature T = t/20 
(200 K) that is high compared with phase transitions @ 
but that is below the pseudogap temperature T*. We 
have verified that they remain essentially unchanged [24| 
as T is decreased towards the range where experiments 
have been performed 0, . 

Our results for 15% doping (Fig. ^a)) are in remark- 
able agreement with ARPES data on reduced e-doped 
cuprates 0, Q when we use U = 5.75*. In ARPES, su- 
perconductivity leads only to a 2 meV shift of the lead- 
ing edge [|| so it would not be observable on the scale 
of Fig. We see that A < (kp,oj) is peaked near (ir,0) 
and (7r/2, 7r/2) and is pseudogaped at hot spots. Smaller 
values of U can be below the critical U c necessary for 
spin fluctuations to build up in the presence of frustrat- 
ing t' and t". At U — 5.25*, for instance, sharp quasi- 
particle peaks appear everywhere along the Fermi sur- 
face (not shown here). For larger U, the spectral weight 
A < (kp,uj) is strongly suppressed near (ir/2,ir/2) because 
scattering by spin fluctuations is stronger and extends 
over a broader region in k space 

M- For 10% e-doped 
cuprates, on the other hand, the best fit to experiments 
is found for U = 6.25i (Fig.^b)), not for U = 5.75t. At 
U = 6.25* the weight A< (hp, u>) near (7r/2, 7t/2) is shifted 
away from the Fermi energy (pseudogap) , in good agree- 
ment with the experimental results 0, while at U = 5.75t 
(not shown here), a pseudogap would not appear near 
(7r/2, 7r/2) 25]. An increase of U with decreasing dop- 
ing is also expected from the fact that U — 6.25* is not 
large enough to produce the observed |26J Mott gap at 
half-filling. Our calculated A < (fc, oS) along some selected 
directions also agrees well with ARPES data 0- Along 
(0,0) — (7r,7r) (Fig. Hie)) a sharp peak approaches the 
Fermi energy near kp, while along (0.657T, 0) — (0.65-7T, w) 
(Fig.QJd)) the peak does not go through the Fermi energy 
and seems to bounce back near kp. The peak position 
where bouncing occurs may be used to obtain a rough es- 
timate of pseudogap size Apq. For 15% and 10% dopings 
we find about 0.1 and 0.3 eV, respectively. 

Calculated Fermi surface plots A < (fc, 0) are obtained 
by integrating A(k,oj)f(uj) in a window (—0.2*, 0.1*), 
comparable to that of experiment @- The overall fea- 
tures are not sensitive to this choice for the parameters 
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FIG. 5 


: (color) Fermi surface plot for (a) U = 


6.25i, 1 



T = t/20 (b) U = 5.75t, n = 1.15, T = i/20 (c) [7 = 5.75t, 
n = 1.20, T = t/50 in the first quadrant of the Brillouin zone. 
See the corresponding experimental plots in Fig. 3 of Ref. 0. 
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used here. For 10% doping (Fig. Ufa)), strong AFM fluc- 
tuations cause the suppression of spectral weight not only 
at hot spots, but also along a large segment of the Fermi 
surface near (ir/2,ir/2). As a consequence, at 10% dop- 
ing the Fermi surface plot looks as if it was composed of 
pockets near (7T,0). As shown in Fig. |3^b) , for 15% dop- 
ing the spectral weight at uj = is strongly suppressed 
by AFM fluctuations only at hot spots. All these results 
are observed in ARPES data 0, for e-doped cuprates. 
These features follow from k dependent scattering rates 
due to AFM fluctuations. 

b. Pseudogap temperature and quantum critical point 
Within TPSC, it has been predicted analytically that the 
pseudogap in the spectral function occurs at T* when 
points of the two-dimensional Fermi surface can be con- 
nected by the AFM wave vector [27j and when £ beg ins to 
exceed £ th = frv F /irk B T at these points [IE El Ha • The 
pseudogap is a precursor of the long-range AFM state in 
two dimensions. Decreasing T below T* increases £ be- 
yond £th and does not produce large changes in the spec- 
tra [24[. In our calculations, the AFM-induced T* ends 
at the same value of doping as the QCP where T = 
two-dimensional antifcrromagnctism disappears. For x 
below the QCP, the T = magnetic order can be com- 
mensurate or incommensurate although at T* the fluc- 
tuations are commensurate. The two-dimensional QCP, 
which may be masked by a superconducting Kosterlitz- 
Thouless transition [2^ or by three-dimensional ordered 
phases, is near 18%, much smaller than the naive perco- 
lation threshold @- At the QCP we find w/T scaling in 
the self-energy and £ ~ with dynamical exponent 

z = 2. This asymptotically implies £ < £th and hence 
no pseudogap at the QCP. Fig. [5f c) illustrates our pre- 
diction that hot spots in A < (k, 0) disappear beyond the 
QCP even at very low temperature. We are also predict- 
ing that the pseudogap observed in ARPES energy distri- 
bution curves at low temperature should persist almost 
unchanged up to the T* indicated by the filled circles in 
Fig. As temperature increases above T* , the pseudo- 
gap in energy distribution curves begins to disappear and 
red segments in Fig. |3 will enlarge gradually until one 
recovers a continuous Fermi surface at around IT* . 

The empty circles in Fig. |3| are the observed pseudo- 
gap temperatures T* extracted from optical conductivity 
data by Onose et al. O- This T* follows closely the T* 
that we predict for ARPES and scales with the observed 



FIG. 3: Pseudogap temperature T* (filled circles denote cal- 
culated T* , empty ones experimental data extracted from op- 
tical conductivity [H^). Empty triangles are experimental 
Neel temperatures Tjv. The samples are reduced Q. 



Tjv curve (approximately two times larger) supporting 
the AFM fluctuation origin of the pseudogap in the nor- 
mal state of e-doped cuprates. Onose et al. 0] found 
that A PG and T* satisfy the relation A PG ss l0k B T*. 
Our estimation of the ratio of Ape to T* leads to 10 
and 7.5 for 10% and 15% dopings, respectively, in rough 
agreement with the above experiment. Note that T* is 
considerably decreased compared with the mean-field Tn 
by both Kanomori screening and thermal fluctuations, 
two effects taken into account by TPSC. 

c. Antijerromagnetic correlation length To demon- 
strate convincingly that our picture is consistent, we 
need to show that our results also agree with the ex- 
perimentally observed AFM £. Mang et al. Q and Mat- 
suda et al. Q have obtained neutron scattering results 
for £(T) at various dopings for oxygenated samples of 
Nd2- a; Ce a ;Cu04. They find an exponential dependence 
of £(T) on inverse temperature. This is illustrated on an 
Arrhenius plot in Fig. |H alon g wi th our calculated £(T). 
Analytical results for TPSC jlfjj indeed show that the 
pseudogap occurs in the renormalized classical regime 
where the characteristic spin fluctuation energy is less 
than the thermal energy and where £(T) ~ exp(C/T) 
with C being generally a weakly temperature dependent 
constant. The temperature at which A < (k 1 uj) has been 
calculated, T — i/20 denoted by an arrow in Fig.0J is still 
in the exponential regime. The value of £(T) at that tem- 
perature for optimal doping is 20 to 30 lattice constants 
(long but finite). For 10% doping (not shown here), the 
calculated £(T) is also in good agreement at high tem- 
peratures but the agreement becomes less satisfactory at 
low temperatures since the quantitative details of TPSC 
become less accurate for very large £. Also, one should 
allow for some small uncertainty in the doping |8( or in 
the parameters of the Hubbard model because neutron 
experiments were performed on oxygenated (as grown) 
samples while ARPES was done on reduced ones |29| . 
We predict that dynamic neutron scattering will find 
overdamped spin fluctuations as we enter the pseudo- 
gap regime and that the characteristic spin fluctuation 
energy will be smaller than k B T whenever a pseudogap 
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FIG. 4: Semi-log plot of the AFM correlation length (in units 
of the lattice constant) against inverse temperature (in units 
of J = 125 meV) for x =0.14, 0.15, and 0.18. Filled symbols 
denote calculated results and empty ones experimental data of 
Refs.0 and ® (x = 0.15). The arrow indicates T = t/20 
where most spectral quantities in Figs. Q and |2] are calculated. 



is present. Equality will occur well above T* . 

d. Long-range order In the extended TPSC ap- 
proach |3j AFM fluctuations can both help and hinder 
d-wave superconductivity, depending on the strength of 
the pseudogap. For the case at hand, a finite T c ap- 
pears below 18% doping, and increases with decreasing 
doping. We cannot study at this point the effect of the 
appearance of long-range Neel order on superconductiv- 
ity but the estimated value of T c near optimal doping is 
roughly the experimental result 0, 0] . By adding hop- 
ping t z = 0.03i in the third direction, we also find a Neel 
temperature that is in fair agreement with the experi- 
mental data by Mang et al. ; .8j for x > 0.14. 



In contrast to e-doped cuprates, in h-doped systems 
a pseudogap opens up near (ir, 0) while spectral weight 
persists near (tt/2,it/2). TPSC can also qualitatively re- 
produce this result but the AFM correlation length in 
this case disagrees with experiment. A different short 
range mechanism can come into play in the strongcou- 
pling limit more appropriate for h-doped cuprates [5j- 

In conclusion, we have shown that ARPES spectra 
0, Q and the antiferromagnetic correlation length ob- 
tained by neutron scattering 0, m e-doped cuprates 
can all be theoretically explained in detail within the 
Hubbard model in the weak to intermediate coupling 
regime. We have made several predictions for future ex- 
periments. The physical picture of the normal state that 
emerges is that of electrons on a planar lattice scatter- 
ing off quasi-static antiferromagnetic fluctuations. These 
have large phase space in two dimensions so that quasi- 
particles in regions of the Fermi surface that can be con- 
nected by the antiferromagnetic wave vector do not sur- 
vive scattering by these low-energy bosons 0, H3 • 
The electron-doped cuprates appear to be the first high- 
temperature superconductors whose normal state lends 
itself to such a detailed theoretical description. 
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